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THE SCATTERING OF OBLIQUELY INCIDENT PLANE WAVES 


FROM A CORRUGATED CONDUCTING SURFACE 
David M. Le Vine 

ABSTRACT 

A physical optics solution is presented for the scattering 
of plane waves from a perfectly conducting corrugated sur- 
face in the case of waves incident from an arbitrary direction 
and for an observer far from the surface. This solution is 
used to compute the radar cross section of the surface in the 
case of baekscatter from irregular (i. e. , stochastic) corru- 
gations and is used to point out a correction to the literature 
on this problem. 

An interesting feature of the solution is the occurrence 
of singularities in the scattered fields. These singularities 
appear to be a manifestation of focussing by the surface at 
its "stationary 1 ' points. Whether or not the singularities 
occur in the solution depends on the manner in which one 
rescripts the analysis to the far field. 
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INTRODUCTION 


The scattering of electromagnetic radiation from irregular surfaces is a 
problem of relevance to a number of remote sensing applications and in particu- 
lar to sensing of the ocean surface. An important special case, which can be 
examined for insight into the m general problem, is that of plane waves inci- 
dent on corrugated surfaces \ surfaces whose ordinate is a function of only one 
dimension). It is the purpose of this paper to present a physical optics solution 
for this special case for an observer in the far field of a conducting surface and 
for plane waves incident obliquely on the surface corrugations. This is a three 
dimensional problem which reduces to the two dimensional form (plane of inci- 
dence perpendicular to the corrugations) as a special case. 

This work was motivated by the physical optics treatment of the scattering 
of plane waves from an irregular surface by Kodis (1966). This solution was 
formulated for the most general case, employing a three dimensional dyadic 
Green's function and an arbitrarily directed plane wave incident on the arbitrary 
surface, z « Z(x, y). The two dimensional case was deduced as a special case; 
unfortunately, an important step seems to have been omitted in making the tran- 
sition to two dimensions and as a result the physical optics solution obtained for 
this special case exhibits an incorrect dependence on frequency and distance. 
These results appear to have propagated into the literature (Barrick, 1968). 

The same error also occurs elsewhere in the literature on scatter from 
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irregular surfaces, appearing in the text by Beckmann and Spizzichino (Bedkmann 
and Spizzichino, 1966). 

The solution to be presented here is for the transitional case of plane 

I 

•.vaves incident obliquely on a corrugated surface. An interesting feature of 
this solution are singularities which can occur depending on the manner in which 
the "far field" approximation is made. If this approximation is employed in a 
minimal way only to justify asymptotic evaluation of integrals, singularities 
appear in the solution for the scattered fields. The singularities disappear in 
the limiting case of an observe infinitely far above the surface, in which case 
the solutions reduce to those obtained by making the far field approximations in 
a conventional manner (i. e. , to simplify the argument of the exponents before 
evaluation of the integrals). The singularities are not peculiar to the case of 
plane waves or to two dimensions, but rather they appear to be a manifestation 
of focussing by the parabolic arcs with which the asymptotic evaluation of integrals 
effectively represent the surface. Roughly speaking, the surface is replaced in 
the asympi-ocic limit by an infinitesimal "mirror" which, with proper curvature, 
can focus the incident radiation at the observation point. This is a characteristic 
of the asymptotic approximation to the integrals and, to the extent that the asymp- 
totic approximation is part of the physical optics approach, the singularities are 
inherent to the analytical technique employed rather than a manifestation of the 
particular problem treated here. 



SOLUTION 


It is assumed that a piano wave, e(r, v) =■ E 0 'e exp (jk • r) is incident on me 
perfectly conducting surface, z = Z(y), where Z(y) is an arbitrary function of y 
and 'e > is a unit vector in the direction of the electric field vector. A physical 
optics solution is to be obtained: that is, Maxwell's equations are "integrated" 
by means of the Kirchoff (tangent plane) approximation and then the integrals are 
evaluated asymptotically in the "high frequency" limit. In order to handle waves 
incident at oblique angles the solution will be formulated initially as a throe di- 
mensional problem. Thus, letting G 0 (r/r') denote the free space dyadic Green's 
function and h(r, v) the magnetic field intensity on the surface, one obtains the 
following form for the (complex representation of) the scattered electric field: ' 


s (r, p) ~ j com (jjf G 0 (r/r) • [n X h(r, i«)l ds' 


( 1 ) 


where 


G„(r/r') 




e jkR 

4aR 


( 2 ) 


and R = |r - r'| is the distance from a point on the surface, r\ to the observer 
at r, and k = 2jtp/c. Employing the Kirchoff approximation, one obtains the fol- 
lowing form for the vector product of h(r, v) with the surface normal, n: 


1i X h(r, p) = * < k(n • c) - e(n • h) 

cum ^ 


)j. e i kk * 


(3) 
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Substituting Equation 2 and 3 into Equation 1, yields the following result fdr the 
scattered fields: 



At this point the results are general in the sense that uo assumptions have 
been made on the direction of k or on the surface, z; and the results are equiva- 
lent to those presented by ICodis (1966). The omission cited above occurred in 
the specialization of these results to corrugated surfaces (z = Z(y)) and to waves 
incident normal to the corrugations; in particular, the x' -coordinate, which in 
tills case no longer appears in the expression for the plane waves or the surface, 
was neglected, and Equation 4 was treated as if the integration was only over y'. 
Unfortunately, in beginning with a three dimensional formulation, one must treat 
the x' -integration explicitly regardless of the plane of incidence. In order to 
avoid tills extra integration, one would have had to formulate the solution initially 
as a two dimensional problem (i. e. , in terms of the two dimensional Green’s 
function). 

It is possible, in the process of making the transition from the three dimen- 
sional case to the strictly two dimensional case, to handle the intermediate prob- 
lem of plane waves obliquely incident on a corrugated surface. This is the case 




to be examined here. Thus, assume a corrugated surface, z « 2(y), and a plane 
wave incident from an arbitrary direction. Consider the x'-integration first and 
make the following transformation of coordinates (Senior, 1959): 

(y'-y ) 2 +(Z(y) -z) 2 = a 2 


x'~x = asinlif?) 


one obtains: 


R = a cosh? 
dx' = a cosh? d? 

and the expression for the scattered field becomes: 


e s (r,r) = -j2kE Q J e j k I k yy' + ^zZ<y')3 T x (y') dy' 


where: 


\ x (y') = 


F(?, y') exp {jka [cosh? + k x sink?] J d? 


Now consider the case in which ka » 1. (a will be large if the observer is f ar 
enough above the mean surface.) Then, neglecting terms In F(?, y ! ) which are 
inversely proportional to lea and representing I x (y') by the first term in its asymp 
totic expansion in ka, one obtains: 

!*</>'= ( 8 ) 

where S' = \A 2 + k 2 a(y ! ) and ? Q = tanh" 1 (~k x ). One will recognize the expression 
in braces in Equation 8 to be the asymptotic form of H^ l \ka). Denoting this 
asymptotic form symbolically by a tilda over the Hankel function, H^(ka^), one 
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obtains upon substituting Equation 8 into Equation G: 


c s (r, v) = -j2kE 0 e^ kk * x J y')[” H^^ka)] c jk M' + M(v')l dy > 


(9) 


Notice that when k x = 0 (i.e. , the plane of incidence is the y-z plane), 

7 0 = 0 and a(y) = a{y). This case is just the two dimensional problem and in 
this case Equation 9 reduces to the result that one would have obtained by begin- 
ning initially with a two dimensional formulation in terms of the Green’s function, 
j// ^(kR). The implications of the asymptotic evaluation of Equation 7 in the 

more general case, k x ^ 0, is that scattering takes place at such points that the 
distance from the scatter points to the observer is \/ (y 1 - y) 2 + (Z(y) - z) 2 * 
cosh( 7 Q ). This implies that the spherical (angular) coordinates of the line from 
observer to scatter point are such that: sin 0 Q cos 0 O = -k x . 


Assuming that k is large, one can now perform the integration over y 1 asymp- 
totically. Doing so, one obtains the following result for the asymptotic value of 
Equation 9: 


cosa(y n ) e 


jMy n ) 




n (r ^ p (y ^ 

s °an? n n 0 


( 10 ) 


where: 


cos 0j + 5 0 cosP 0 


= k y y n + K zcy n ) + 5 0 «(y n ) + k x x 

6(y n ) 

W 


cos a(y n ) cos 2 (a + j3 0 ) 

[1 - (VR)(VR)] • [fc(n • c) - c(n -1c)] 
6 


a la) 
(lib) 
(11c) 
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Notice that a(y n )/8 0 is just the distance from the observer to the n-th scatter 
point. When this distance is large in comparison to R c (y n )/e(y n ), Equation 30 
indicates that the field scattered from each stationary point decreases as the 
squareroot of the distance from the scatter point to the observer as is charac- 
teristic of scattering from a two dimensional object. The field measured by 
the observer as predicted by Equation 10 appears to come from many separate 
scattering centers, one for each y n , and the magnitude of the radiation scattered 
from each such point to the observer depends on the squareroot of the distance 
from the scattering point to the observer, a(y n )/5 0 , on the squareroot of the 
radius of curvature, R c (y n ), at the scatter point and on the relative orientation 
of the incident ray, the observer, and the slope of the surface at the scattering 
point (i.e. , on 0j, J3 0 anda(y n )). 

The integration over the x 1 - and y' -coordinates have resulted in a pair of 
restrictions on the coordinates of the scatter points possible for a given incident 
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wave and fixed observer. The x’ -integration leads to the restriction, sin 0 6 • 
cos tj> 0 = -k x , where 0 o and 0 O are the spherical coordinates (measured at the 
observer) of the line from observer to the scatter point. It follows that for each 
0 Q there are two possible configurations: <p o - ±cos" 1 (-k x /sin 0 Q ). These two 
possibilities are illustrated in Figures 1 and 2. They correspond to the case of 
"forward" scatter (k y in the reflected ray has the samo sign as in the incident 
ray) as shown in Figure 1 and the case of "back" scatter as illustrated in Figure 
2. In the case of back scatter, reflection takes place at points on the surface at 
which a(y) > 0 whereas forward scatter occurs at points for which o!{y) ^ 0. When 
k x = 0, then ~ ±7r/2 and in this case both incident and reflected waves are per- 
pendicular to the surface corrugations. This is the degenerate case correspond- 
ing to the two dimensional problem. There also are a forward and back scatter 
possibility in this case. (See Figure 3. ) In forward scatter the reflected and 
incident waves both propagate in the same direction with respect to the y-axis, 
and in the case of back scatter they propagate in opposite directions with respect 
to the y-axis. 


The condition imposed by the y f -integration is that — - = 0 = k v + k, tan a(y) 

3y * 

3a 


4- \/k x + k^ — . This condition is equivalent to the requirement that the local 
angle of incidence and reflection are equal. In order to see that this is so, it 

rw 

is convenient to define angles 0 O and to be the (acute) angles that the projec- 
tions of the reflected and incidentrays onto the y-»- plane make With the vertical 
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(i.e. , the z-axis), Then, in terms of , Fj anti a, the condition 
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comes sin (fy, - a) - sin (a - fy). This implies that tho angle of incidence equals 
the ang’e of reflection as measured in the plane of the projection. Since tho plane 
of incidence also contains the local normal (and the reflected i*ay) it is also true 
that the anglo of incidence equals the angle of reflection as measured in tho piano 
of Incidence. That is, the points at which reflection takes place are "specular 11 
points. (Tho fact that the normal to the surface, the incident ray and tho reflected 
ray all lie in tho same plane— the plane of incidence— is a consequence of tho 
boundary conditions imposed by the tangent plane approximation. ) 


The amplitude of the radiation scattered from each stationary point and 
measured at the observer is given by Equation 10, An interesting feature of 
this amplitude is the singularity which occurs at <=(y n ) a(y n )/5 Q = ~R c (y,, ). Since 
e(.V„> < 0 for physically meaningful geometry (i. e. , waves incident toward the 
surface so that lCj, < 0), it follows that R c (y n ) < 0 at scatter points associated 
with the singularity. Thus, at these points the surface is concave toward the 
observer. It appears that the singularities are the manifestation of focussing 
which can take place at the surface due to the combination of tangent plane ap- 
proximation (which preserves the phase structure of the incident plane wave) 
and the asymptotic evaluation of the integrals (which effectively represent the 
surface by an arc with radius of curvature, R c (y n )). Since the surface is per- 
fectly conducting, the result is that each scatter point behaves as if a plane wave 
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were incident on a concave or convox mirror, depending on tho sign of the tad.ku? 

of curvature. When the mirror is concave, focussing can occur. For example, 

consider the case of two dimensions, k x = 0. In this case, 5 0 = 1 and 0 O - 0 Q and 

= 0j whore fi Q and aro the obtuse angles associated with P 0 and p ( , respcc- 

tively. In this caso, the requirement that « 0 yields tho relationship, 0 o = 

oy 

2a - 0[, and it follows that e(y n ) = 2sec(a - 0 O ). Consequently, the singularity 
occurs whenever a(y n ) = 2sec(a - 0 o )/R c (y n ), or in terms of the acute angle, 0 Q , 
at a(y n ) = - 2 sec(a - (T 0 )/R c (y n ). Now a(y n ) is just the distance along the scattered 
ray £r..»n$ the scatter point to the observer; and if one imagines a concave spherical 
mirror to be located at the scatter point with radius of curvature, -R c (y n ), and 
an axis which coincides with the normal to the surface, then -2 sec(a - 0 o )/R c (y n ) 
is just the distance along a reflected ray from the mirror to the focal plane. (See 
Jenkins and White, 1957.) 

The results shown in Equation 10 were obtained by making a somewhat modi- 
fied "far field’* approximation. In particular, the factor, R = |r - r’l , wnskeptin 
an arbitra.,/ rorm in the exponential, exp <jkR), but was treated as a large num- 
ber with comensui'ate simplifications being made whenever this was expeditious 
only in multiplicative factors. This is in contrast to the more conventional ap- 
plication of the far field approximation in which R is approximated by the constant 
and linear terms in its binomial expansion for use in the exponential prior to eval- 
uation of the integrals. (For example, see Kodis, 1966.) If one follows this latter 
approach, no singularities occur and the form obtained is just Equation 10 in the 
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limiting case o£ e(y„)a(y n )/5 0 » R c (y n ) for all R c (y n ). That is: 


E _ / costly ; ^ 

F n (y n ) cos a(y n ) / — 

all y n v cos °i ’ h 5 o cos P<) V a (y n y 5 o 


cos a(y ) 

— r t==== VKJyT) (12) 


Apparently the focussing is overlooked in making the far field approximation in 
tilts conventional form because this effectively removes the observer to infinity, 
beyond all possible focal points. It is interesting to note that the singularities 
also occur in the treatment of finite sources (Le Vine, 1974) and therefore are 
not peculiarities of plane waves. 


RADAR GROSS SECTION OF A STOCHASTIC SURFACE 

As a final point, a comparison will be made of the radar cross section in two 
dimensions of a stochastic surface based on the results derived here (Equation 
12) and those based on the computation of Kodis (Kodis, 1966). In keeping with 
the treatment of Kodis, it is assumed that the receiving antenna is sufficiently 
narrow and properly oriented so that scatter only from points for which the incident 
and scattered ray are colinear need be considered. Furthermore, it is assumed 
that the k0(y n ) are uniformly distributed over 2yr and that they are independent 
of R c (y n ), (This assumption amounts to assuming that the incident radiation 
scatters incoherently from the surface.) With these assumptions, the radar 


cross section, o° , for a unit length of surface and based on Equation 12 is: 

I / Lim fl .Efl\ 
a 0 = -( 2t rp > 

L \> -+ oo L E o J / 

= irN<R c (y n )> 


(13a) 

03b) 








where L is the length of surface illuminated, the pointed brackets denote an en- 
semble average, N is the number of scatter points per unit length of the illuihi- 
nated portion of the surface, and p = a(y)/5 0 is the distance from the observer 
to the surface. Equation 13a is just the two dimensional adaptation of the usual 
formula for radar cross section (Kerr, 1951} Skolnik, 1970), In comparing 
Equations 12 and 13 with those in Kodis (Kodis, 1966, Equations 14 and 26) 
notice that there is no frequency dependence in the results derived here, whereas 
Kodis' results are proportional to k. Notice also that power, in the two dimen- 
sional case as presented here, decays as 1/p (i. e, , as a cylindrical wave) as op- 
posed to inversely as the square of distance as occurs in the three dimensional 
case and in Kodis' results (Kodis, 1966, Equations 14 and 15). The three dimen- 
sional character of Kodis' results are the consequence of the omission discussed 
above. 
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Figure 3. Scatter in the Two Dimensional Case: a) Forward Scatter; 

b) Back Scatter 
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